In-silico enhanced animal study of pulmonary artery pressure sensors: assessing hemodynamics using computational fluid dynamics.

To assess whether in-silico models can be used to predict the risk of thrombus formation in pulmonary artery pressure sensors (PAPS), a chronic animal study using pigs was conducted. Computed tomography (CT) data was acquired before and immediately after implantation, as well as one and three months after the implantation. Devices were implanted into 10 pigs, each one in the left and right pulmonary artery (PA), to reduce the required number of animal experiments. The implantation procedure aimed at facilitating optimal and non-optimal positioning of the devices to increase chances of thrombus formation. Eight devices were positioned non-optimally. Three devices were positioned in the main PA instead of the left and right PA. Pre-interventional PA geometries were reconstructed from the respective CT images, and the devices were virtually implanted at the exact sites and orientations indicated by the follow-up CT after one month. Transient intra-arterial hemodynamics were calculated using computational fluid dynamics. Volume flow rates were modelled specifically matching the animals body weights. Wall shear stresses (WSS) and oscillatory shear indices (OSI) before and after device implantation were compared. Simulations revealed no relevant changes in any investigated hemodynamic parameters due to device implantation. Even in cases, where devices were implanted in a non-optimal manner, no marked differences in hemodynamic parameters compared to devices implanted in an optimal position were found. Before implantation time and surface-averaged WSS was 2.35±0.47 Pa, whereas OSI was 0.08±0.17, respectively. Areas affected by low WSS magnitudes were 2.5±2.7 cm2, whereas the areas affected by high OSI were 18.1±6.3 cm2. After device implantation, WSS and OSI were 2.45±0.49 Pa and 0.08±0.16, respectively. Surface areas affected by low WSS and high OSI were 2.9±2.7 cm2, and 18.4±6.1 cm2, respectively. This in-silico study indicates that no clinically relevant differences in intra-arterial hemodynamics are occurring after device implantation, even at non-optimal positioning of the sensor. Simultaneously, no embolic events were observed, suggesting that the risk for thrombus formation after device implantation is low and independent of the sensor position.

Comparison of hemodynamics in biological surgical aortic valve replacement and transcatheter aortic valve implantation: An in-silico study.

OBJECTIVES: In aortic valve replacement (AVR), the treatment strategy as well as the model and size of the implanted prosthesis have a major impact on the postoperative hemodynamics and thus on the clinical outcome. Preinterventional prediction of the hemodynamics could support the treatment decision. Therefore, we performed paired virtual treatment with transcatheter AVR (TAVI) and biological surgical AVR (SAVR) and compared hemodynamic outcomes using numerical simulations. METHODS: 10 patients with severe aortic stenosis (AS) undergoing TAVI were virtually treated with both biological SAVR and TAVI to compare post-interventional hemodynamics using numerical simulations of peak-systolic flow. Virtual treatment procedure was done using an in-house developed tool based on position-based dynamics methodology, which was applied to the patient's anatomy including LVOT, aortic root and aorta. Geometries were automatically segmented from dynamic CT-scans and patient-specific flow rates were calculated by volumetric analysis of the left ventricle. Hemodynamics were assessed using the STAR CCM+ software by solving the RANS equations. RESULTS: Virtual treatment with TAVI resulted in realistic hemodynamics comparable to echocardiographic measurements (median difference in transvalvular pressure gradient [TPG]: -0.33 mm Hg). Virtual TAVI and SAVR showed similar hemodynamic functions with a mean TPG with standard deviation of 8.45 ± 4.60 mm Hg in TAVI and 6.66 ± 3.79 mm Hg in SAVR (p = 0.03) while max. Wall shear stress being 12.6 ± 4.59 vs. 10.2 ± 4.42 Pa (p = 0.001). CONCLUSIONS: Using the presented method for virtual treatment of AS, we were able to reliably predict post-interventional hemodynamics. TAVI and SAVR show similar hemodynamics in a pairwise comparison.

CT-Based Analysis of Left Ventricular Hemodynamics Using Statistical Shape Modeling and Computational Fluid Dynamics.

Background: Cardiac computed tomography (CCT) based computational fluid dynamics (CFD) allows to assess intracardiac flow features, which are hypothesized as an early predictor for heart diseases and may support treatment decisions. However, the understanding of intracardiac flow is challenging due to high variability in heart shapes and contractility. Using statistical shape modeling (SSM) in combination with CFD facilitates an intracardiac flow analysis. The aim of this study is to prove the usability of a new approach to describe various cohorts. Materials and Methods: CCT data of 125 patients (mean age: 60.6 ± 10.0 years, 16.8% woman) were used to generate SSMs representing aneurysmatic and non-aneurysmatic left ventricles (LVs). Using SSMs, seven group-averaged LV shapes and contraction fields were generated: four representing patients with and without aneurysms and with mild or severe mitral regurgitation (MR), and three distinguishing aneurysmatic patients with true, intermediate aneurysms, and globally hypokinetic LVs. End-diastolic LV volumes of the groups varied between 258 and 347 ml, whereas ejection fractions varied between 21 and 26%. MR degrees varied from 1.0 to 2.5. Prescribed motion CFD was used to simulate intracardiac flow, which was analyzed regarding large-scale flow features, kinetic energy, washout, and pressure gradients. Results: SSMs of aneurysmatic and non-aneurysmatic LVs were generated. Differences in shapes and contractility were found in the first three shape modes. Ninety percent of the cumulative shape variance is described with approximately 30 modes. A comparison of hemodynamics between all groups found shape-, contractility- and MR-dependent differences. Disturbed blood washout in the apex region was found in the aneurysmatic cases. With increasing MR, the diastolic jet becomes less coherent, whereas energy dissipation increases by decreasing kinetic energy. The poorest blood washout was found for the globally hypokinetic group, whereas the weakest blood washout in the apex region was found for the true aneurysm group. Conclusion: The proposed CCT-based analysis of hemodynamics combining CFD with SSM seems promising to facilitate the analysis of intracardiac flow, thus increasing the value of CCT for diagnostic and treatment decisions. With further enhancement of the computational approach, the methodology has the potential to be embedded in clinical routine workflows and support clinicians.

CT-Based Simulation of Left Ventricular Hemodynamics: A Pilot Study in Mitral Regurgitation and Left Ventricle Aneurysm Patients.

Background: Cardiac CT (CCT) is well suited for a detailed analysis of heart structures due to its high spatial resolution, but in contrast to MRI and echocardiography, CCT does not allow an assessment of intracardiac flow. Computational fluid dynamics (CFD) can complement this shortcoming. It enables the computation of hemodynamics at a high spatio-temporal resolution based on medical images. The aim of this proposed study is to establish a CCT-based CFD methodology for the analysis of left ventricle (LV) hemodynamics and to assess the usability of the computational framework for clinical practice. Materials and Methods: The methodology is demonstrated by means of four cases selected from a cohort of 125 multiphase CCT examinations of heart failure patients. These cases represent subcohorts of patients with and without LV aneurysm and with severe and no mitral regurgitation (MR). All selected LVs are dilated and characterized by a reduced ejection fraction (EF). End-diastolic and end-systolic image data was used to reconstruct LV geometries with 2D valves as well as the ventricular movement. The intraventricular hemodynamics were computed with a prescribed-motion CFD approach and evaluated in terms of large-scale flow patterns, energetic behavior, and intraventricular washout. Results: In the MR patients, a disrupted E-wave jet, a fragmentary diastolic vortex formation and an increased specific energy dissipation in systole are observed. In all cases, regions with an impaired washout are visible. The results furthermore indicate that considering several cycles might provide a more detailed view of the washout process. The pre-processing times and computational expenses are in reach of clinical feasibility. Conclusion: The proposed CCT-based CFD method allows to compute patient-specific intraventricular hemodynamics and thus complements the informative value of CCT. The method can be applied to any CCT data of common quality and represents a fair balance between model accuracy and overall expenses. With further model enhancements, the computational framework has the potential to be embedded in clinical routine workflows, to support clinical decision making and treatment planning.